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ABSTRACT 

The statistics of large-separation gravitational lensing are a powerful tool to probe 
mass distributions on the scale of galaxy clusters. In this paper we refine the analysis 
of our survey for large-separation (> 5") lensed FIRST quasars (Ofek et al. 2001) 
by estimating the magnification bias and the source redshift distribution. Finding 
no large separation lens among 8000 likely quasars in that sample, implies an upper 
bound on the lensed fraction of 3.7 x 10 -4 at 95% CL. From a published deep 1.4 GHz 
radio survey of the Hubble Deep Field, and corresponding optical searches for faint 
quasars, we calculate a lower limit to the 'double flux magnification bias' affecting 
our radio-optically selected sample, of B > 1.1. From the four-colour information in 
the SDSS Early Data Release, we calculate the photometric redshift distribution of a 
sample of FIRST quasar candidates and compare it with the redshift distribution from 
the FIRST Bright Quasar Survey. We find that the median redshift of the quasars in 
our sample is about 1.4. With these new results, we find that for all plausible cos- 
mologies, the absence of lensed quasars in our survey is consistent with a model based 
on an empirical, non-evolving, cluster mass function, where clusters are represented 
by singular isothermal spheres. On the other hand, comparison of our results to the 
lensing predictions of published iV-body-ray-tracing simulations (Wambsganss et al. 
1995, 1998) rejects the COBE normalised fi = 1 CDM model at 99.9% confidence. 

Key words: cosmology: gravitational lensing - galaxies: clusters: general - quasars: 
general 



1 INTRODUCTION 

The statistics of gravitational lensing can provide a powerful 
probe of the geometry and the mass content of the universe 
out to large redshifts (e.g. Refsdal 1964; Press & Gunn 1973; 
Turner, Ostriker, & Gott 1984). Lensed quasar statistics in 
the angular image separation ranges expected from galaxy- 
mass lenses have been probed by several surveys (e.g., Maoz 
et al. 1993). However, systematic searches for lensed quasars 
with multiple image separations of 5" A8 <i 30" expected 
from groups or clusters of galaxies have been carried out only 
recently (e.g., Maoz et al. 1997; Marlow et al. 1998; Phillips 
et al. 2001a; Phillips, Browne, & Wilkinson 2001; Ofek et 
al. 2001; Keeton & Madau 2001; Phillips et al. 2001b) 

The incidence of large-separation gravitational lensing 
depends on several factors, including the cosmology, the 
present-day cluster mass function and its redshift evolution, 



* e-mail: eran@wise.tau.ac.il 



and most importantly, the mass profile and the substruc- 
ture within galaxy clusters (e.g., Bartelmann, Steinmetz, & 
Weiss 1995). There are several approaches for predicting the 
optical depth to large-separation lensing in a given model 
scenario: 

(i) Ray tracing through iV-body simulations (e.g., Cen 
et al. 1994; Wambsganss et al. 1995; Wambsganss, Cen, 
& Ostriker 1998). This accounts for the complicated large- 
scale, multi-plane lens structure, and for the substructure in 
clusters of galaxies. Wambsganss et al. (1995; 1998) predict 
an optical depth for lensing with AG > 5" of about 3 x 10 -3 , 
for a z ~ 3 source in an Qo = 1 Cold Dark Matter (CDM) 
cosmology. The problem of substructure, in particular, was 
analysed by Bartelmann & Weiss (1994), and Bartelmann 
et al. (1995). They found that numerically modeled clus- 
ters with substructure are about two orders of magnitude 
more efficient in producing highly magnified images (e.g., 
long arcs) than spherically symmetric clusters. 

(ii) Taking the observed cluster mass function and 
adopting a smooth and symmetric mass profile to calculate 
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the optical depth (e.g., Maoz et al. 1997; Ofek et al. 2001). 
Maoz et al. (1997) used the NFW profile (Navarro, Frenk, 
& White 1996, 1997) to show that the cluster lensing opti- 
cal depth can be in the range t(A6 > 5") ^ 10~ 5 to 10~ 2 , 
where r depends mostly on the adopted mass-concentration. 
In Ofek et al. 2001 (hereafter Paper I), we estimated the op- 
tical depth for singular isothermal sphere (SIS) lenses using 
the non-evolving Girardi et al. (1998) cluster mass function 
and found t{A6 > 5") w 3 x 10~ 4 , for an Q = 0.3, Ov = 0.7 
cosmology and a source redshift of z 30urce — 1. 

(iii) Modifying approach (ii) by assuming the cluster 
mass function is given by the Press-Schechter formalism 
(Press & Schechter 1974), that by itself depends on cosmol- 
ogy. For example, Li & Ostriker (2002) have followed this ap- 
proach, and assumed a generalized NFW density profile and 
z source = 1- Their predictions for the optical depth span val- 
ues in an even larger range, 7 x 10~ 7 ^ r(A8 > 5") ^ 10~ 3 , 
depending in part on the details of the adopted mass pro- 
file. The impact of the cosmological model on the cross sec- 
tions may seem surprisingly large, about one or two orders 
of magnitudes. This dependence enters through the differ- 
ing redshift evolutions of the cluster mass function. Other 
examples of the Press-Schechter approach to cluster lensing 
statistics are Narayan & White (1988), Flores & Primack 
(1996), Wyithe, Turner, & Spergel (2001), and Sarbu, Rusin, 
& Ma (2001). 

In light of these wide-ranging predictions, empirical 
constraints seem in order. Without any confirmed large- 
separation quasar lenses, we can only place an upper limit 
on the cross section which, in turn, can reject part of the 
mass-profile and cosmology parameter space. 

In Paper I we described a survey for large separation 
gravitational lensing among quasars drawn from the FIRST 
radio catalog. In that survey, we selected radio sources from 
the FIRST catalogue (1999, July 21 version; Becker, White, 
& Helfand 1994; White et al. 1997) with optical counter- 
parts in the APM catalog (McMahon, & Irwin 1992), with 
positional coincidence better than 2. "5 in both the O and E 
bands, with colour index O — E < 2 mag, and that are point- 
like on at least one of the O or E POSS-I/UKST plates, as 
determined by the APM classification algorithm. There are 
12, 576 quasar candidates passing the above criteria. Based 
on the initial results from the FIRST Bright Quasar Survey 
(FBQS; Gregg et al. 1996; White et al. 2000), we estimated 
in Paper I that there are about 9100 quasars in our sample 
out of the 12, 576 candidates. 

However, with the more recently published results of the 
FBQS (Becker et al. 2001), which extends the spectroscopic 
follow-up from 17.8 to 19 mag, we can revise the fraction 
of quasars among the sources fainter than E — 17.8 mag in 
our sample from 90% to 80%, or about 8000 quasars instead 
of 9100 quasars. Specifically, Equation 1 in Paper I can be 
revised to: 

_ J -0.9389 + 0.046£ + 0.0027£; 2 , E < 18 . . 
tqso ~ \ 0.8, E > 18 ( ' 

Follow-up observations of all pairs of radio-optical 
quasar candidates with 5" < A9 < 30" showed that none 
of the 15 candidate pairs are lensed quasars (for details see 
Paper I). For 8000 quasars in the sample, this implies an 
upper limit of 3.7 x 10~ 4 (95% Poisson statistics confidence 
level; CL) on the lensing fraction in this survey. 



However, in order to constrain the mass function and 
mass profile of clusters from these results, some further prop- 
erties of our survey need to be known: (i) the double flux 
magnification bias of our survey, resulting from the radio 
and optical flux limits involved in the sample selection. This 
magnification bias reflects the over- (or under) representa- 
tion of lensed objects in our survey, compared to purely geo- 
metrical cross-section calculations; and (ii) the redshift dis- 
tribution of quasars (i.e., the potential sources) in our sam- 
ple, which determines the effective pathlength of our lens 
search. 

The purpose of the present paper is to improve those 
two important constraints and provide a more reliable up- 
per limit on the observed optical depth to large separation 
lensing. 



2 THE SURVEY DOUBLE FLUX 
MAGNIFICATION BIAS 

As described in Paper I, the magnification bias depends 
strongly on the faint-end number counts (in the radio and 
optical) of quasars with similar selection criteria to those 
used in our survey. However, there are presently no large- 
area very deep radio-optical surveys. The best available data 
set is the combination of radio and optical surveys in the 
Hubble Deep Field North (HDF-N) region: a deep 1.4 GHz 
radio survey conducted with the VLA by Richards (2000). 
This survey has a flux limit (5a) of 40 ^iJy in the cen- 
tral 8', degrading to about 300 /iJy at 18' from the field 
centre. In the optical, a 1 deg 2 region around the HDF- 
N was searched for quasars down to B = 22 mag using 
multi-colour selection (Liu et al. 1999; Vanden Berk et al. 
2000). Liu et al. required somewhat more restrictive color 
criteria than ours: B — R < 0.8, or U — B < —0.4 and 
B - R < 1.1. Impey & Petry (2001) matched the Richards 
(2000) radio sources with the optically detected quasars in 
this field, and identified three radio-optical quasars within 
the colour range, 0.46 < B — R < 0.75, that match our 
survey criteria. With sufficient lensing magnification these 
quasars would have been detected by our survey. These data 
on the HDF-N constitute faint-end counts that can be com- 
bined with our FIRST- APM candidates, to find the double 
flux magnification bias. 

The HDF implies radio counts of 190tioo radio-optical 
sources per deg 2 per mjy, at 0.12 mjy in the 1.4 GHz-band, 
and optical counts of 3. 51^9 radio-optical sources per deg 2 
per mag at O = 20.3 magnitude. The new radio number- 
count point includes three radio-optical quasars with optical 
counterpart brighter than B = 22, and the optical point is 
based on one faint (fainter than the APM 100% complete- 
ness limit of 19.75; Caretta et al. 2000) radio-optical quasar 
with radio counterpart brighter than the Richards (2000) 
survey limit, as given above. Figure 1 shows the differential 
number counts in the 1.4 GHz-band for the Paper I sample 
quasars (circles), and for the deep radio-optical search for 
quasars described above (single square). The dashed line is 
the interpolated number-counts in the 0.1-1 mjy range. Fig- 
ure 2 is the same as Figure 1, but for the optical O-band. 

Although it is reasonable to assume that there is a large 
number of radio- faint (fainter than 0.1 mjy) quasars consti- 
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Figure 1. Differential number counts in the 1.4 GHz-band for 
our sample quasars (circles) , and one square representing the deep 
radio-optical search for quasars described above. The dashed line 
is the interpolated number-counts in the 0.1-1 mjy range. 
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Figure 2. Same as Figure 1, but for the O-band 



tuting the radio-quiet quasar population, we only have di- 
rect observational constraints on the number counts down 
to 0.1 mjy (Figures 1-2). Therefore, our estimate for the 
double-flux magnification bias will only be a lower limit. 
Furthermore, due to the small area of the HDF-N deep radio 
observations, the number counts at the intermediate flux lev- 
els (between FIRST and HDF-N) are not measured directly, 
and we simply interpolate the number count distribution 
over these gaps. 

Drawing on this new information we have calculated the 
double flux magnification bias in the FIRST-APM sample 
through a Monte-Carlo simulation, as follows. 
We constructed a list of radio-optical quasars that mimics 
the observed number count distribution of the O-magnitudc 
and of the 1.4 GHz flux shown in Figures 1-2. For quasars 
above the FIRST-APM flux limits we drew at random ob- 
jects from our list of quasar candidates, corrected for the 
fraction F qso of these objects that are quasars, and for the 



incompleteness of the APM survey (Caretta et al. 2000), by 
replicating objects in the magnitude-range 19.5-21.5, until 
we had a list of about 8000 radio-optical sources. At the faint 
end we added objects with radio flux and optical magnitude 
drawing from a flux probability distribution obtained by in- 
terpolating the number counts in Figure 1. For each of these 
objects we then drew an O-magnitude from the distribu- 
tion given in Figure 2. Note that the last step assumes that 
there is no correlation between the optical and radio prop- 
erties of quasars, which is true at least for objects brighter 
than 1 mjy and O ~ 21 (Paper I). The number of randomly 
selected objects in this step was chosen by scaling the num- 
ber of HDF-N faint sources by the ratio of the total area of 
the FIRST-APM survey to the area of the HDF-N faint-end 
number-count survey described above. 

With the list described above, we are now in a position 
to estimate the double-flux magnification bias, assuming a 
SIS lens mass profile. Our FIRST-APM survey will detect 
a lensed quasar as such only if both lensed images (bright 
and faint) are brighter than the radio 1 mjy flux limit, and 
are brighter than the O-magnitude limit of 21.2, which is 
the weighted average magnitude limit of the APM (Caretta 
ct al. 2000). The magnification/de-magnification of the faint 
image, A-, in a SIS lens is given by: f3 = 1/(1+j4_), where /3 
is the impact parameter in units of the Einstein ring radius. 
By comparing the magnitude and flux of each object in the 
Monte-Carlo list with the survey's flux limits, we calculated 
for each object, i, the minimum magnification needed (or 



maximum de-magnification allowed), Ai 



and the cor- 



responding fli^max, needed for the detection of an image in 
both the radio and optical bands in our survey (i.e., the 
smaller among the optical and radio values of $ 

,max WaS 

chosen) . 

The cross-section for strong lensing is defined as the 
area with (3 < 1, and the probability for a given impact pa- 
rameter (3 scales like the area, /3 2 . To estimate the double- 
flux magnification bias, B, we summed (3l, max for all the 
objects in our Monte-Carlo list and divided the results by 
the number of quasars in our sample, i.e., 8000. The re- 
sulting ratio reflects the over- (or under) representation of 
lensed objects in our survey, compared to a purely geomet- 
rical cross-section calculation. 

Using the above method, we obtain a lower limit of 
B ~Z 1.1. This limit changes to 0.9 if we take into account 
the la Poisson statistics (lower confidence interval of the 
number counts), and to 1.3 if we assume that the optical 
searches for quasars described above are 50% incomplete 
(see Impey & Petry 2001) due to the different colour selec- 
tion process used by Liu et al. (1999) and Vanden Berk et al. 
(2000). The B > 0.9 limit we give above assumes that the 
number-counts are almost flat from O ~ 19 to O ~ 20, and 
then go to zero. The actual value of B could be larger due to 
the probable existence of a yet-fainter radio-optical quasar 
population. Hence, we consider B > 1.1 to be a conserva- 
tive limit. For comparison, in Paper I we roughly estimated 
B by assuming that the product of the optical and radio 
number-flux relations is a power-law with a cutoff at a flux 
that is 0.5 — 4 orders of magnitude lower than our survey's 
flux/magnitude limits. Based on this, we estimated B to 
be in the range of 1 — 40. Since the HDF-N data extend the 
depth of our optical source counts by about 1 mag, the lower 
limit on B found in this paper is near the low-end estimate 
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of Paper I, but is now on a firmer observational basis. The 
present method we have used to calculate the bias is supe- 
rior to the Borgeest, Linde, & Refsdal (1991) approach we 
used in Paper I, since it takes into account the actual shapes 
of the separate optical and radio number count distributions 
instead of replacing them by single power-law distributions. 



3 REDSHIFT DISTRIBUTION OF THE 
FIRST/ APM SOURCES 

From the first major public release of the Sloan Digital Sky 
Survey - Early Data Release (SDSS-EDR; York et al. 2000; 
Stoughton et al. 2002), we constructed a sample of 443 point- 
like SDSS objects that: were selected as quasar-candidates 
in the SDSS-EDR based on their colors; are found within 
l."2 from a FIRST radio source; and are brighter than 
g' — 21 mag. Their magnitudes were corrected for Galac- 
tic extinction (Schlegel, Finkbeiner, & Davis 1998). All but 
two sources in this list have g' — r' < 2 mag, so they are good 
representatives of our survey. Based on the SDSS quasar se- 
lection criteria (e.g., Newberg & Yanny 1997; York et al. 
2000; Schneider et al. 2002), and on the FBQS (Gregg et 
al. 1996; Helfand et al. 1998; Becker et al. 1998; White et 
al. 2000; Becker et al. 2001), we know that the fraction of 
quasars in this sample is J; 75%. 

From the quasar redshift- colour relation, given by 
Richards et al. (2001a, Table 3), we estimated the photo- 
metric redshifts for our 443-object sample, by a simple \ 2 
minimization (Richards et al. 2001b) between the predicted 
and the observed SDSS colours. The photometric-redshift 
distribution of these objects is shown as bars in Figure 3. 
The dashed line shows the redshift distribution of quasars 
in the FBQS faint extension (E < 19 mag; Becker et al. 
2001). The photometric-redshift distribution is not changed 
considerably if only FIRST sources with SDSS counterparts 
brighter than g' = 19 mag are taken. However, our photo- 
metric redshift estimate of quasars has a considerable scat- 
ter: 30% of the redshifts have errors, Az > 0.25. To estimate 
the effect of this scatter on the photometric redshift distri- 
bution, we have used the 898 quasars with known redshifts 
and SDSS colours given by Richards et al. (2001a) to cal- 
culate the distribution of differences between our redshift 
prediction and the observed redshifts in each 0.25 redshift 
bin. We convolved this scatter distribution as a function of 
redshift with the FBQS faint-extension redshift distribution 
(E < 19 mag; Becker et al. 2001) and the result is shown 
as the dotted line in Figure 3. Ideally, we would expect that 
the photometric redshift distribution (bars) be similar to 
the redshift distribution of FBQS quasars convolved with 
the photometric-redshift scatter (dotted line). However, our 
photometric-redshift distribution has an excess of quasars at 
z ~ 3, relative to the FBQS faint-extension redshift distri- 
bution The median redshift for the sample of FBQS faint- 
extension quasars is z = 1.2, and the median redshift we get 
from the photometric redshift down to limiting magnitude 
g' — 21 is z — 1.4. There are two possibilities to explain the 
discrepancy: (i) only a fraction (<; 75%) of the sources in our 
sample are real quasars. Therefore, the non-quasar objects 
could contaminate the redshift distribution in some system- 
atic way; (ii) the photometric-redshift procedure of Richards 
et al. (2001b) was developed based mostly on radio-quiet 



0.18r 

0.16 
0.14 
0.12 
0.1 



S0.08 

20.06 
Q_ 

0.04 



0.02 



Figure 3. The redshift distribution of FIRST quasar candidates 
down to g 1 = 21 mag (bars). The redshift distribution was cal- 
culated using a simple colour-redshift relation. The superposed 
vertical-lines on the bars, are the Poisson error-bars due to the 
finite sample size. The dashed line is the redshift distribution of 
FBQS faint-cxtension (E < 19 mag) quasars. The dotted line is 
the same distribution, but convolved with the photometric red- 
shift scatter, as a function of redshift. 



quasars, but our survey contains mainly radio-loud quasars. 
It is known that there are some optical spectral differences 
between radio-loud and radio-quiet quasars (e.g., Brother- 
ton et al. 2001). Such differences could bias the Richards 
et al. (2001b) photometric-redshift procedure. Therefore we 
adopt the FBQS faint-extension redshift distribution for our 
calculations. This is a conservative choice, as the lower me- 
dian redshift implies a smaller path length for lensing. 



4 DISCUSSION 

In the following sections we use these two new pieces of infor- 
mation, B > 1.1 and the redshift distributions, to estimate 
the optical-depth to lensing in a simple SIS model, and to 
compare our results with Qo = 1 CDM TV-body ray-tracing 
simulations. 



4.1 A SIS model 

To get some simple estimates of the lensing probability 
in our survey, we calculated the cross section as a func- 
tion of redshift, assuming SIS lenses and the Girardi et al. 
(1998) mass function of groups and clusters with no evo- 
lution. We found in Paper I that the effective dimension- 
less density of lenses (as defined by Turner et al. 1984) is, 

F = 167r 3 n (-^) 3 (^f) 4 = 0.0156±g;°^ 7 , where c is the 
speed of light, Ho the Hubble constant, an the line-of-sight 
velocity dispersion, and no is the comoving number-density 
of lenses. The errors in F were obtained from Monte-Carlo 
simulations using the Girardi et al. (1998) mass function pa- 
rameters and their errors (for details see Paper I). In Paper I 
we assumed a single representative redshift of z = 1 for all 
the quasars. With our new results presented above, we re- 
calculate the total optical-depth for lensing, by integrating 
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the optical-depth over redshift, with the weighting function 
given by the FBQS faint-extension redshift probability den- 
sity distribution (dashed-line in Figure 3). The optical depth 
r z , to a source at redshift z, is given by: 

f Z ev-i , ^3 fDiD ls \ 2 1 cdt 

where Di, Di s , D s are the angular-diameter distances for the 
observer-lens, lens-source and observer-source, respectively, 
zi is the lens redshift, and Ro = c/Hq. The quantity cdt/dzi 
is calculated in the Friedmann-Lemaitre-Robertson- Walker 
geometry (e.g., Fukugita et al. 1992). The total weighted 
optical depth is given by: 




where P(z) is the normalized weighting function given by 
the FBQS quasar redshift distribution. We find 

r 4.6±2:o x 1(r4 > Q o = 0.3; Ov = 0.7 
r = I 2A±\ 8 x 10~ 4 , n = 0.3; Ov = 0.0 (4) 
1 1.7±o;f x 10~ 4 , n = 1.0; n A = 0.0, 



where the errors are the la errors due to the uncertainty 
in F. This is ~ 1.7 times greater than estimated in Pa- 
per I, where z = 1 was assumed for all sources. Using the 
photometric-redshift distribution from Figure 3 (bars) in- 
stead of the FBQS faint- extension redshift distribution, the 
predicted r would further increase by a factor of about 
1.6. The difference is due to the excess of quasars with 
z ~ 3 in our photometric-redshift distribution, combined 
with the fact that the optical-depth increases with red- 
shift. We conservatively adopt the optical-depth in Equa- 
tion 4, as a lower limit on r. These optical-depths need 
to be multiplied by the completeness factor that reflects 
the range of separations probed. As shown in Paper I, this 
completeness factor is about 0.66 for a SIS model. Finally, 
the lower limit on the magnification (B > 1.1) determines 
the expected probability for lensed quasars in our sample: 
Piens > 4.6±2'.o x 10~ 4 x 0.66 x 1.1 = 3.3±l B B x 10~ 4 (for 
fio = 0.3, Q A = 0.7). 

The observational upper bound on the lensed fraction 
of 3.7 x 10~ 4 (95% CL) from our survey (Paper I, and small 
revisions in this paper), is therefore consistent with this 
simple model (Piens)- Since the (Oo = 0.3, SIa = 0.0) and 
(fio = 1.0, fl A = 0.0) models predict a lower optical depths, 
they too are consistent with the data. Within this simple 
model, for an flo — 0.3, Qa = 0.7 cosmology and a non- 
evolving Girardi et al. (1998) mass function, we can in turn 
set an upper limit to F < 0.018 with a 95% CL. Models 
incorporating evolution of the mass function will reduce the 
optical-depth, therefore increasing the upper limit on F. 

4.2 Comparison with ray-tracing simulations 

Wambsganss et al. (1995, 1998) performed ray tracing for 
< 2 < 3 through a three-dimensional mass distribution 
obtained from an flo — 1 COBE normalized CDM iV-body 
simulation. They found that the optical depth for strong 
lensing (i.e., for forming multiple images) with image sepa- 
rations > 5" of sources at z ~ 3 is 3x 10 . Our observations 
can test this model. Figure 2 in Wambsganss et al. (1995) de- 
scribes the splitting probability for image pairs with separa- 
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tions greater than 5" and magnitude differences smaller than 
1.5, as a function of the source redshift. To obtain the ex- 
pected optical-depth for our survey, we integrate their split- 
ting probability with a weight function given by the FBQS 
faint-extension source redshift distribution presented in Fig- 
ure 3. We obtain r = 1.14x 10~ 3 (or r = 1.5x 10~ 3 using the 
photometric-redshift distribution) . This optical depth needs 
to be multiplied by the completeness factor for the range 
of separations probed and the magnification bias of our sur- 
vey. Wambsganss et al. (1995) also show (their Figure 3) the 
strong-lensing probability as a function of the image separa- 
tion for various source redshifts. From this figure we estimate 
the incompleteness due to the separation range probed in 
our survey, P(A > 30")/P(A > 5") ~ 0.25 (This separation 
incompleteness is actually slightly redshift dependent, and 
changes from ~ 30% at z s — 1 to ~ 20% at z s — 3). Taking 
together r, the separation incompleteness and the magnifica- 
tion bias we have found here, implies a probability for lensing 
in our survey P C dm ~ 1.14 x 10~ 3 x 0.75 x 1.1 = 9.4 x 10~ 4 
based on the Wambsganss et al. simulations. 

We note that the magnification bias we found in this 
paper is based on a SIS flux ratio distribution - an assump- 
tion that is not necessarily compatible with the Wambsganss 
et al. (1995, 1998) simulations. However, Figure 5 in Wamb- 
sganss et al. (1998) shows that the magnification probabil- 
ity found in their simulations is not very different from the 
P(A) ~ A expected from the SIS model we assumed in 
our magnification-bias calculation. Based on the calculation 
described above, 7.5 lenses are expected in our survey for an 
Q.q = 1 CDM model. We can therefore reject the Wambs- 
ganss et al. (1995, 1998) model at 99.9% CL. We note that, 
since Figure 2 in Wambsganss et al. (1995) gives the splitting 
angle probability for image flux ratios smaller than 1.5 mag, 
our conclusion is conservative. Note, that the simulations of 
Wambsganss et al. use the COBE normalization, as = 1.05, 
which predicts an excess of high-mass clusters at the present 
epoch, and slower evolution of the mass function. In this 
sense, our results do not reject the fi = 1 model per se, but 
the combination of the cosmology and the high normaliza- 
tion. 

The results presented here constitute one of the tight- 
est and well-defined limits on large separation lensing to 
date. Clearly, we have now reached a level where the 
non-detections are becoming interesting cosmological con- 
straints. It is probable that actual cases of large-separation 
gravitationally lensed quasars will be found soon using one 
of the large area surveys (e.g., SDSS, 2DF; see Croom et 
al. 1998). However, in order to use such large separation 
lenses to constrain the cluster mass function and mass pro- 
file parameter space, we need to: (i) understand the selection 
effects and effectiveness of these surveys for large separation 
lensing; (ii) have a realistic estimate of the cross section for 
large separation lensing, taking into account the substruc- 
ture of clusters. These points will be addressed in future 
papers. 

To summarize our main results, we have derived a lower 
limit to the magnification bias in our survey, B > 1.1, and 
we have also found that the median photometric-redshift 
of quasars in our sample is z source ~ 1.4. Using the red- 
shift distribution from the FBQS faint extension, we deter- 
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mine that F < 0.018 with 95% CL (assuming a SIS model 
and an fio = 0.3, Qa = 0.7 cosmology). Our non-detection 
of lensed FIRST quasars is consistent with expectations, if 
clusters can be represented by a non-evolving population of 
SIS masses with the mass function of the observed cluster 
population. If so, moderately larger surveys will discover the 
first examples of large separation quasar lensing. Our survey 
already has the ability to reject some models with concrete 
predictions, namely the flo = 1 COBE-normalized CDM 
model, whose excess of power on large scales is well known. 
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